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Abstract. We report on searches for B° — > fj,~ , B d — > decays and b — > s/z + /x _ 

transitions in exclusive decays of B mesons. Using 2 fb" 1 of data collected by the CDF II 
detector we find upper limits on the branching fractions B(B® — > /i + /i~) < 5.8 x 10 -8 and 
3{B ( l -> < 1-8 x 10" 8 at 95% confidence level. Using 924 pb _1 of data we measure the 

branching fractions B(B+ -> ^+^~K + ) = (0.60 ± 0.15 ± 0.04) x 10~ 6 , B{B° d -> ^+ K*°) = 
(0.82 ±0.31 ±0.10) x 10~ 6 and the limit B(B° -> fi + ^ (f>)/B(B° -> J/ij)<t>) < 2.61(2.30) x 10 -3 
at 95(90)% confidence level. 



1. Introduction 

In the standard model (SM) both the decay of a B® d meson into two muons and the decay of a 
b quark into an s quark and two muons are highly suppressed. They require a Flavor Changing 
Neutral Current (FCNC) process which can only occur through higher order Feynman diagrams. 
Processes beyond the standard model can contribute to those decays and significantly modify the 
decay rates. Any sizable deviation of the measurement from the SM prediction would indicate 
new physics phenomena. In these proceedings we present the current searches of the CDF 
Collaboration for the branching ratios of the rare decays B® s « — > [1] and B — > fi + [i~h 

[2], where B stands for B + , B d or B®, and h stands for K + , K*° or respectively. 

2. Search for the rare decays B^ s d ^ — > 

The SM predictions for the branching ratios of the decays B® d — > are of the order of 

1(T 10 - 10 -9 [3]. These low rates can not be probed at the Tevatron yet. However, new physics 
models like supersymmetry can enhance those ratios to a level that can be observed by the 
Tevatron experiments. 

The analysis uses 2 fb _1 of data [4], selected by the CDF dimuon trigger. The CDF muon 
system consists of several different components, most notably the central muon system (CMU), 
covering the pseudorapidity range \r]\ < 0.6, and the central muon extension (CMX), covering 
0.6 < \t]\ < 1.0. The dimuon trigger system requires at least one muon in the CMU. The 
two resulting trigger possibilities, CMU-CMU and CMU-CMX, have different acceptances and 
detector efficiencies, so that they will be treated separately during the analysis. 
After the application of some baseline cuts, the selection of B® d — > candidates uses 

vertex displacement, transverse momentum and isolation variables that are combined in a neural 
network (NN) to form a single variable, discriminating between signal and background. For 
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Figure 1. The invariant dimuon mass distribution versus the NN output for events satisfying 
the baseline cuts in the CMU-CMU mode (left). The NN output bin NN> 0.995 is the most 
sensitive search area. Two candidates are found in the B® mass range, five candidates in the B d 
mass range. The right plot shows the dimuon mass distribution for events with NN> 0.995. 



network training B® d — > signal events from simulation and background events from data 

sidebands are used. The invariant dimuon mass distribution is shown in Fig. 1. 

The optimal search strategy is determined using the a priori expected 95% C.L. upper limit 
as figure of merit. The search window is divided into 15 regions - 5 mass bins times 3 NN 
output bins. The observed number of events is compatible with the estimated background. 
Combinatoric background is estimated by sideband interpolation and the contribution from 
misidentified hadronic B® d decays is obtained by a combination of data and simulation. The 
limit on the branching fraction is obtained by normalizing to the well-measured B + — ► J/ipK + 
decay. After combining all 15 regions, the CMU-CMU and CMU-CMX modes and using 
the world average value for B(B + — ► J/ r ipK + ) [5], we obtain the 90%(95%) C.L. limits of 
B{B° S -> m+aO < 4 -7 x 10~ 8 (5.8 x 10~ 8 ) and B{B Q d -> < 1.5 x 1(T 8 (1.8 x 10" 8 ). These 

results are the most stringent limits to date. 

3. Search for the rare decays B + — > fi + fi~~K + , B° — » fi + fi ~K*° and B° — > ~<j) 
The transition b — * s[i + fi~ is sensitive to new physics since physics beyond the SM can manifest 
itself in significant deviations from SM transition rate predictions. The analysis uses a data 
sample corresponding to 924 pb _1 . Each branching fraction is determined by the corresponding 
resonant B — > J/i^h mode for normalization, resulting in a cancellation of many systematic 
uncertainties. After some baseline cuts and the removal of resonant dimuon contributions the 
selection is optimized using the following variables: the proper decay time significance ct/o~ c t, 
the pointing angle from the B meson candidate to the primary vertex and the isolation of the 
B meson candidate. The figure of merit in the optimization process is S/y/S + B, with S the 
expected signal yield, determined from the normalization modes, and B the expected background 
yield. The following background sources are considered: misidentification of hadrons as muons, 
reflections between the three modes, and combinatorial background. 

The three final invariant mass distributions are shown in Fig. 2. Each mode shows an 
excess in the signal mass region. The significance of the excess is determined by calculating 
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Figure 2. Invariant mass spectra for the rare modes. The black curve illustrates the expected 
shape after background estimation and determination of the signal yield. The dotted part of the 
curve indicates that it is not expected to describe the low mass background, containing partially 
reconstructed B decays. 



the probability for the expected background to fluctuate into at least the number of observed 
events. From the event counts, together with the normalization sample, the relative branching 
fractions between signal mode and normalization mode are obtained. The absolute branching 
fractions are calculated with the world average branching fraction values for the normalization 
modes [5]. Table 1 summarizes the results. 



Table 1. Summary of the B — * /U + /i h results. Listed for each decay mode are the number of 
observed signal events, the number of expected background events and the signal significance. 



Mode 


B+ -> n + ii-K + 


B° d -> n+^K*° 


B° s -> tj,+fi-<f> 


N observed events 


90 


35 


11 


N expected BG events 


45.3 ±5.8 


16.5 ±3.6 


3.5 ±1.5 


Gaussian significance 


4.5 


2.9 


2.4 


Rel B ± stat ± syst x lO" 3 
Abs B ± stat ± syst x 1(T 6 
Rel B 90%CL limit xKT 3 
Rel B 95%CL limit xl(T 3 


0.59 ±0.15 ±0.03 
0.60 ±0.15 ±0.04 


0.62 ±0.23 ±0.07 
0.82 ±0.31 ±0.10 


1.24 ±0.60 ±0.15 
1.16 ±0.56 ±0.42 

2.30 

2.61 



The results for B{B + — > /x + pT K + ) and B(B® — > fi + K*°) are in good agreement with 
previous measurements [6, 7] and SM predictions [8]. The limit on B(B® — > n + fi~(fi) is the most 
stringent to date. 
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